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We analyze the trajectories for multifield DBI inflation, which can arise in brane inflation models, 
and show that the trajectories are the same as in typical slow roll inflation. We calculate the power 
spectrum and find that the higher derivative terms of the DBI action lead to a suppression of the 
contribution from the isocurvature perturbations. We also calculate the bispectrum generated by 
the isocurvature perturbation, and find that it leads to distinctive features. 



INTRODUCTION 

Scalar field theories with non-canonical kinetic terms 
provide novel realizations of the inflationary paradigm 
|EJ. One interesting class of such models which have 
been studied extensively in recent years is DBI inflation 
0, |j| , characterized by kinetic terms which arise from the 
Dirac-Born-Infeld (DBI) action. A particularly appeal- 
ing phenomenological feature of DBI inflation is that it 
can lead to strong and unique non-Gaussian signatures 
in the Cosmic Microwave Background (CMB) |,[§. The 
DBI action is ubiquitous in string theory, e.g., as an effec- 
tive theory for worldvolume degrees of freedom on branes. 
When embedded in brane inflation fjj, the inflaton field 
in DBI models can be given a natural geometrical inter- 
pretation as the position of a D-brane in extra dimen- 
sions. DBI inflation arises when the D-brane moves in 
a highly warped region of the internal space where the 
speed limit is small, and reduces to the usual slow roll 
brane inflation (with canonical kinetic term) when the 
brane is moving non-relativisitically with respect to the 
local warp factor. Since the position of the brane in each 
compact direction is described by a scalar field, brane 
inflation is naturally a multifield inflationary model. 

Multifield models are characterized by their trajecto- 
ries in field space and can in general be decomposed into 
an adiabatic field, which parameterizes motion along the 
trajectory, and isocurvature fields, which describe the 
directions perpendicular to the trajectory. Features in 
the trajectory, such as a sharp turn, can convert isocur- 
vature perturbations into adiabatic/curvature perturba- 
tions (even on superhorizon scales) and can give rise to 
interesting features in the primordial power spectrum and 
non-Gaussianity. 

In this paper we study the effects of multiple fields in 
DBI inflation. In particular, we study the multifield DBI 
trajectories and show that they are identical to the usual 
slow-roll multifield case in which the trajectory is domi- 
nated by the field with the largest slope of the potential. 
We calculate the power spectrum for multiple DBI fields 
in the limit the trajectory makes a sharp turn and show 
that the contribution of the isocurvature perturbations 



to the power spectrum is suppressed by the sound speed. 
Finally, we calculate the non-Gaussianity in the sharp 
turn limit and find that the non-Gaussianity is domi- 
nated at leading order in the sound speed by the usual 
single field DBI contribution, but has new multifield fea- 
tures at sub-leading order. We conclude by commenting 
on the impact of our results for DBI model building. 



MULTIFIELD EQUATIONS OF MOTION 

Consider a 10-dimensional warped throat background 
with the metric 

ds 2 = f- 1/2 (y)g^dx»dx» + f 1 / 2 {y)~g mn dy m dy n (1) 

common to type IIB string compactifications ^, Q , where 
f{y) is called the warp factor of the throat which can in 
principle depend on all of the coordinates of the internal 
space y. 

A D3-brane in this background is described by (to low- 
est order in string coupling and to all orders in a'), 
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where the warp factor is rescaled f{<j>i) = T^ 3 f(y((/)i)), 
and the real canonical scalar fields associated with the 
motion of the brane are given by 
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where i = 1, ...6. The potential V(4n) can arise, for ex- 
ample, from interactions with £>3-branes, D7-branes, or 
from the breaking of the local isometries of the compact 
space, and we will leave it to be unspecified for the mo- 
ment. 

We will define the sound speed during inflation to be 
the inverse of the "Lorcntz factor" of the DBI action for 



spatially homogeneous fields, 



(4) 



We will be interested in the small sound speed limit 
c s « 1 where the non-Gaussianity is observable. In the 
following we will use the convention that the yi measure 
the distance from IR "tip" of the throat and that the D3 
brane moves towards the tip, so </>j < 0. 

Consider a FRW universe with four-dimensional met- 
ric, 



ds 2 = -dt 2 + a(t) 2 ^2 dx i- 
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The Friedman equation and equations of motion for the 
fields fa are (with H = -) 
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Distributing the time derivative, the equation of motion 
(0) can also be written 
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Clearly the equation of motion for a homogeneous scalar 
field with a canonical kinetic term is obtained from 
(0),(|) in the limit c s -> 1. 

We define the multifield DBI inflationary parameters 

as a 
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When these parameters are much smaller than one, the 
equations of motion ^) take a form similar to that of a 
slowly rolling field 

Y,, taxiOijrjij + e 
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+ drb, (y(^)+/ _1 (c S -l)) 

« 3H4h + c s d, k (V{4>i) + r\c s - 1)) = , (11) 

where tan^y » d^V/d^V to leading order in the DBI 
inflationary parameters (p[jlO|). 

As a specific example, for the "standard case" of an 
AdS warp factor that depends only on one of the fields 



FIG. 1: The trajectory of a multiple field inflationary system 
can be decomposed into an "adiabatic" field a with compo- 
nents along the trajectory and an "entropy" field s orthogonal 
to the trajectory. 



f{<j>i) = A/0 4 and a separable potential of mass terms 
V = |X) TO f^fj ^ can be shown that for small sound 
speed c s « 1 and a large mass hierarchy (m, ^> nij for 
some the multifield DBI inflationary parameters are 
of order 0(e). 



TRAJECTORIES IN MULTIFIELD BRANE 
INFLATION 

To simplify our analysis we will restrict ourselves to 
a two field model ((pi,<p2), but it is straightforward to 
generalize our analysis to any number of fields. 

We will parameterize the classical trajectory by an 
"adiabatic field" a that represents the component of the 
field motion along the trajectory ||] (see Figure |l|), 



(cos 9)4>i + (sin( 
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The angle 9 parameterizes the angle the classical trajec- 
tory makes with one of the field directions (here chosen 
to be (pi) and should not be confused with the angular 
position of the D3 brane in the compact space, which 
is parameterized by the fields <pi. The "entropy field" s 
transverse to the classical trajectory gives rise to isocur- 
vature fluctuations which are given by, 



ds = — (sin + (cos9)S(f>2 
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By definition, the entropy field is constant, s = 0. 

Using this parameterization, we can rewrite the exact 
equations of motion (@) as (in particular, the multifield 
DBI parameters are not necessarily assumed small for 
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this expression), 
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It is clear from (g_5() that the angle of the trajectory has 
a "fixed point" trajectory in field space 9^0 defined by 
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Furthermore this fixed point is stable to leading order 
when the slope of the potential is positive and dominates 
over the slope of the warp factor since small variations 
89 are driven to zero (notice that a < 0). This fixed 
point can also be seen from the equations of motion ( p"l| ) 
in the "DBI slow roll regime" (e.g. when the multifield 
DBI parameters (§10) are small), 
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We see then that being in the DBI slow roll regime is 
equivalent to being at the stable "fixed point" of the tra- 
jectory, thus DBI slow roll is an attractor solution. 

Let us examine ([l6]) in more detail. First, we note 
that for an inflationary solution we require that the po- 
tential energy dominates over the kinetic energy, e.g. 
V{<j>i) > l/(c s /), so that for small c s , V((pi) > c s /f 
is automatically true. The trajectories ( |l6| ) are now the 
same as in standard multifield inflation, in particular, the 
angle of the trajectory is controlled by the ratio of the 
curvatures of the potential: the trajectory follows the di- 
rection with the largest curvature. When the field with 
the largest curvature reaches its minimum the trajectory 
makes a sharp turn in field space, with the sharpness of 
the turn given by the ratio of the curvatures. 

If we express the warped geometry as a cone over an 
angular base space X 5 , 
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we can identify one of the fields <p\ as the radial coor- 
dinate and the other field <p2 as one of the angular co- 
ordinates on the base X 5 of the D-brane. In this case 
we can roughly classify the trajectories based on which 
field dominates at early times, as shown in Figure |[ a 
radially dominated trajectory is when the slope of the po- 
tential in the angular direction <p2 is much smaller than 
the radial direction <p\ and so the trajectory is dominated 
at early times by motion in the radial direction; a diag- 
onal trajectory is when the slopes of potentials of the 
two fields are approximately the same and so the tra- 
jectory is approximately a diagonal line composed of a 
linear combination of the radial and angular directions; 
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FIG. 2: Multifield models have a number of different trajec- 
tories, depending on the curvature of the potential for the 
fields. The sharpness of the turns is controlled by the ratio of 
the curvatures. 



and a angularly dominated trajectory is when the slope of 
the potential in the radial direction (pi is the smallest so 
the trajectory is dominated at early times by motion in 
the angular <p2 direction. Clearly the diagonal trajectory 
does not have a significant turn in field space, and so will 
not contribute to a generation of curvature perturbations 
as discussed above. In fact, the diagonal-type trajectory 
is just a linear combination of the fields and hence can 
be completely described by a single field, and so we will 
not consider this possibility. 



THE POWER SPECTRUM OF MULTIFIELD DBI 
INFLATION 

The equation of motion for the curvature perturbation 
7Z sa C in a multifield model with fields (pi can be solved 
exactly using the SN formalism |), 10 , which states that 
the curvature perturbation is equal to the difference be- 
tween the number of e-folds of the classical trajectory N e 
and the perturbation to the classical trajectory Af e , 



C = dN e = M e - N e = V NiiSfaY 
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where N ti = §|f 



is the derivative of the number of e- 
foldings with respect to the field evaluated when the 
mode exits the horizon <p*. For two fields, the power 
spectrum is then given by, 

4tt 2 P c - (CC> = (^Vi) 2 (<Mi) 

+2{N 1 )(N a )(<p 1 <p 2 ) + {N 2 ?{<t>2<t>2) , (20) 

where we have allowed for cross correlation between the 
fields - for a canonical kinetic term, the cross coupling is 
zero. We see, then, that the power spectrum receives ex- 
tra contributions in multifield inflation, both from cross 
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couplings and from the extra two point correlation func- 
tions of the additional fields. 

We would like to analyze these extra contributions in 
more detail. Using the approximation that the DBI mul- 
tifield parameters ( |l0| ) are small and that the potential 
is separable V(0i,<fe) = ^i(0i) + ^2(^2) and dominates 
the energy density, the number of e-folds from the time 
of horizon crossing to the end of inflation is pi 

1 f e V 1 
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As in jlOj , we can write the curvature perturbation (|19|) in 
a much simpler way by using a different set of multifield 
DBI parameters, 
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where a prime denotes a derivative with respect to the ar- 
gument. Note that e = —H/H 2 = e\ + e 2 . The curvature 
perturbation can now be written, 




with 
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where the superscript 'e' denotes evaluation at the end of 
inflation. In the limit that the trajectory is in the radial 
or angular direction as shown above then at the end of 
inflation the one of the field has settled into its minima 
(say, 4>2 for concreteness) so Z e = V 2 = const. For mass- 
term dominated potentials (or equivalently when the vac- 
uum energy from (j>2 is small V 2 <C V 2 ) then we can take 
Z e = and the expression for the curvature perturbation 
simplifies to include only the values of the potential and 
the slow roll parameters evaluated at horizon crossing, 
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This simple expression for the curvature perturbation will 
be useful later in evaluating the power spectrum. 

The inflationary perturbation 



analytic control over our expressions we will assume that 
the trajectory is highly radial, e.g. 
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tan 9 = — = -fj < c. 

Under this assumption, the adiabatic perturbation is sim- 
ply the perturbation in the <pi direction and perturba- 
tions in the <p 2 direction are just isocurvature perturba 
tions, 
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Expanding the kinetic part of the Lagrangian in ^) to 
quadratic order in 8a and Ss, we find 
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Note that the isocurvature fluctuations (e.g. the fluctu- 
ations in the angular direction) scale differently with the 
sound speed; we will see soon that this has important 
consequences for the two point functions. 

The quantization of the perturbations proceeds as 
usual with 

5<j(t, x) = 



where 
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for the adiabatic perturbation and 
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for the isocurvature modes. The creation and annihi- 
lation operators satisfy the usual commutation relation 
[a(k), a t( k ')] = [6(k),fet( k ')] = (27r) 3 <5 3 (k-k'), «(r,k) 
and v(t, k) are the solutions of the quadratic Lagrangian 
whose normalizations are fixed by the Wronskian condi- 
tions 1 , and r = — \j is the conformal time. 

aH 

It is now straightforward to calculate the two point 
functions, 

(5a(k 1 )«5<7(k 2 )) = (2^ 3 (k 1+ k 2 )^, (31) 



(5«(ki)*a(k 2 )) = (2^) 3 ,5 3 (k 1 +k2) 
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In order to calculate the power spectrum we need to 
evaluate the two point correlation functions for the per- 
turbations of the scalar fields in (pOf). In order to have 



This can be checked by computing the commutator 
[8<j(t, xi),pj CT (r, x 2 )] = i<5 3 (xi - x 2 ), where p Str = is 
the canonical momentum. 
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Here Sa(k) — Sa(r, k)| T _»Oi Ss(k) — Ss(t, k)| T ^ j an d 
the Hubble parameter and sound speed are evaluated at 
the time of horizon crossing. Here we see that the isocur- 
vature fluctuations generated by Ss are suppressed by a 
factor of c s <C 1 compared to the adiabatic perturbations 
generated by 5a. This suppression can be traced back to 
the different c s dependence found in the quadratic La- 
grangian (pcf). We note that the exact same calculation 
follows through if the trajectory is dominated by the an- 
gular direction after the replacement 5<j)\ <-*■ 5<f> 2 . Finally, 
notice also that in the limit of a straight-line trajectory 
there is no cross correlation between the fields, so the 
cross term in the power spectrum (^o|) vanishes. Un- 
fortunately, non-straight line trajectories are beyond our 
analytic control so it is not clear if the cross coupling 
will be significant, although we expect that in the diag- 
onal limit the two point function should reduce to that 
of a simple single field model where the results are well 
known. 

Using the two point functions (31 3^) and the expres- 



sion for the curvature perturbation (25fj, the power spec- 
trum (p0|) becomes, 
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where in the last line we assume that the second term is 
small compared to the first term in the small c s limit. 

Since the contribution of the angular modes to the 
power spectrum is highly suppressed by the sound speed 
we find that multifield DBI reduces essentially to the sin- 
gle field case, in contrast to multiple field slow roll infla- 
tion where additional fields may become important when 
the trajectory makes sharp turns in field space [fto|. 



Multifield Non-Gaussianity 

In the following we study the non-Gaussianities in 
more details, and as it turns out there are some poten- 
tially observable differences from the single field DBI in- 
flation. To compute the non-Gaussianities, we expand 
the DBI Lagrangian to higher order. The leading order 
and subleading order cubic terms are 



r[Sa - a- 2 c 2 s 5a(y5a) 



2c 5 s a 
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The leading contribution to the non-Gaussianity comes 
from the first two terms, and their size is well known 
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(35) 



Since the angular mode 5s is suppressed by a factor of 
yj~c~ s comparing to the radial mode 5a, we see that the 
third term in ([34]) contributes a non-Gaussianity of order 

fNL ~ - (36) 

Cs 

Although it is sub-leading comparing to the first two 
terms, it is still potentially observable by future experi- 
ments if the sound speed c s is small enough. In particu- 
lar, this effect can be larger than the sub-leading effect of 
order \ computed in Q] when the sound speed satisfies 
c s > e. ' 

The calculation of the three-point function is standard, 
see e.g. M for details, 
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here "cyclic." means two other terms by cyclically per- 
muting k\,k 2 ,k 3 . We can see the non-Gaussianity van- 
ishes in the squeezed limit when one of the momentum 



0, since 



du* (V,k3 



— fc 3 and there is factor of ki ■ k 2 . 

This is the same as in single field DBI inflation. 

Going away from the squeezed limit, we can compute 
the above three point function assuming ki , k 2 , k 3 are of 
the same order of magnitude. The non-Gaussianity has 
a very interesting shape as the following 

fc|(fc2 + 3fcifc 2 + A:|)(-ki -k 2 ) 



A{ki,k 2 ,k 3 ) 



(h + k 2 f 



(38) 



-cyclic. 



We plot the non-Gaussianity as A(k\, fc 2 , l)/(fcifc 2 ) fol- 
lowing the convention of 0] in Figure |3| 

We can see from Figure |^ that the non-Gaussianity 
has a very interesting new feature, namely, the sign of 
non-Gaussianity is different in the middle of the "folded 
triangle" limit where fc 3 = k\ + fc 2 from most of the other 
region of the configuration space. For example, one can 
directly check for the configuration of a folded triangle 
k\ = fc 2 = the shape is negative A{k\, k 2 , fc 3 ) — 
— ?i . This feature is not present in any other known 
inflationary models, so it can be used as a distinctive 
signature of multifield DBI inflation. The sign of non- 
Gaussianity determines the sense of skewness in the CMB 
temperature and matter density, and thus this change in 
sign can lead to potentially interesting observational ef- 
fects 2 . However, we should caution that since this is a 



The convention for the sign of /jvx (which characterizes the level 
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FIG. 3: The shape of non-Gaussianity in multifield DBI in- 
flation is shown through a plot of A(ki, ki, l)/(fcifc2) in 
Notice that in the folded triangle limit k\ = ki = kzjl the 
bispectrum is negative, and constitutes a distinctive signature 
of multifield DBI inflation in the small c s limit. The presence 
of opposite signs of the non-Gaussianity may give rise to in- 
teresting observational effects. 




FIG. 4: The negative part of the non-gaussianity in the folded 
triangle limit is shown. 



sub-leading effect, it might be difficult to disentangle the 
signature from experimental data. It would be interest- 
ing to extend this analysis to more than two fields, which 
may lead to additional enhancements or suppressions of 
the multifield non-Gaussianities (l2l. 



Discussion 

We have shown that the two point function of the extra 
"angular" scalar field direction during inflation is sup- 
pressed by a factor of c s compared to the usual single 
field contribution, thus the multifield DBI observables 
simply reduce to the single field case (in the limit where 
one of the fields dominates the trajectory). There have 
been a lot of recent attempts in building single field DBI 
inflation models consistent with all known compactifica- 
tion constraints and precision cosmological observations 
|l3| . In particular, these works find that combining ob- 
servational constraints of the amplitudes and tilt of the 
scalar and tensor perturbations and the primordial non- 
gaussianity together with limitations on the field range 
coming from compactification puts severe constraints on 
the viable parameter space of single field DBI models, 
ruling out the most simple models. 

In general, since we see that the small sound speed in 
multifield DBI inflation suppresses the multifield effects, 
multifield DBI inflationary observables can be well ap- 
proximated by their single field values. Note that this 
also implies that the dramatic 0(1) effects expected at 
the end of multifield DBI inflation due to the inhomoge- 
nous surface of tachyon condensation examined in [fl4| are 
now instead suppressed by 0(c s ) and are subdominant. 

Variations of the basic DBI model, such as its IR ver- 
sion | fl5[ |l6|| , models involving wrapped branes |p"7f or 
different warped geometries, may be able to evade the 
strong constraints of fl3| ; since we have not made ex- 
plicit use of any particular model we expect our results 
to hold in the small sound speed limit of these models as 
well. 

Note: As this manuscript was being prepared, a 
preprint [Ti\ ] appeared which contains some overlap with 
this work. 
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of non-Gaussianity) has been a source of confusion in the litera- 
ture. See JlT| for a discussion. 
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